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Abstract 

We present three preliminary results from SLD on B decays: an inclusive search for the 
process b—>sgluon, a measurement of the branching ratio for the process B -^DDX, and 
measurements of the charged and neutral B lifetimes. All three measurements make use of 
the excellent vertexing efficiency and resolution of the CCD Vertex Detector and the first two 
make use of the excellent particle identification capability of the Cherenkov Ring Imaging 
Detector. The b^sg analysis searches for an enhancement of high momentum charged kaons 
produced in B decays. Within the context of a simple, Jetset-inspired model of sg, a limit 
of B(b-^ sg)< 7.6% is obtained. The B[B — > DDX) analysis reconstructs two secondary 
vertices and uses identified charged kaons to determine which of these came from charm 
decays. The result of the analysis is B(B -> DDX)= (16.2 ± 1.9 ± 4.2)%. The results of the 
lifetime analysis are: r B + = 1.686 ± 0.025 ± 0.042 ps, r B o = 1.589 ± 0.026 ± 0.055 ps and 
t b +/t b0 = 1.061 ±° f 2 l ±0.027. 
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1 Introduction 



Detailed studies of B-hadron decays can provide important tests of the Standard Model. In 
fact, it has been suggested that the existence of several "5-decay puzzles" may be pointing 
the way to physics beyond the Standard Model [Q. The most serious of these puzzles is 
the the low measured value compared to theoretical expectations of the the B semi-leptonic 
branching ratio: @ 

jo ^semi—leptonic /-. \ 

&SL = ^ — = — = (1) 

1 semi—leptonic T J- hadronic T -L leptonic 

Theoretical expectations of Bsl typically have a lower limit of 12.5% ||. As of Summer '97, 
however, the experimental measurements were much lower: the world average of measure- 
ments done at the T(45) was 0.1018 ± 0.0040 and at the Z° was 0.1095 ± 0.0032 §. 

Since r semi _; eptonic is well understood theoretically and Ti eptonic is very small, efforts to 
explain the discrepancy have focused on Thadronic, which can be broken into three parts 
(assuming that Cabibbo suppressed rates are small): 

^hadronic ~ Y(b^cud) + T{b-^ccs) + T(b-^sg) (2) 

Reducing Bsl to the experimentally measured value would require enhancing one or more 
of these components significantly above the expected value. In this paper, we describe 
measurements that bear on each of these three components. 



1.1 SLD Capabilities and Data Set 

The SLD experiment collects Z° decay data from e + e~ collisions at the SLAC Linear Collider 
with a center of mass energy of 91.28 GeV. A full description of the SLD detector may be 
found in [|J. SLD is well-suited for doing precision measurements of inclusive -B-decays 
due to several unique characteristics. The SLC interaction point is small and stable and 
its position is known with an uncertainty of 5 /mi transverse to the beam direction. This 
precise IP is complemented by SLD's high precision CCD vertex detector, VXD3 ||. For 
high momentum tracks the impact parameter resolution is a(r<f)) = 11/iin and a(rz) = 22/xm. 
Multiple scattering adds a momentum-dependent contribution of 33/im/(j?sin 3 / 2 (#)), where 
p is the momentum expressed in GeV/c and 9 is the track polar angle. Note that the 
above describes the performance of the upgraded vertex detector (VXD3) installed prior 
to the start of the 1996 run. For the performance of the vertex detector used before 1996 
(VXD2), see Reference [|J. Another important capability of SLD is the excellent particle 
identification provided by the Cherenkov Ring Imaging Detector |J. K s in the barrel 
region with momentum between 1 and 20 GeV/c are identified with an efficiency of 50% and 
a 7r ± misidentification probability of 2%. Data used in these analyses was taken between 
1993 and 1998. However, not all data has yet been used in all analyses. Table [I] lists the 
different running periods used for each analysis. 
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Running Period Number of Z's Vertex Detector 



b^sg B(B — > DDX) r^+,r R o 



1993-95 150K VXD 2 

1996-98 250K VXD3 

Spring 98 150K VXD3 



Yes No Yes 

Yes Yes Yes 

Yes Yes No 



Table 1: The three running periods used in the analyses are listed along with the vertex 
detector in use for that period and whether the period was used in each of the three analyses 

1.2 Analysis Techniques 

All three of the analyses described in this paper make use of an inclusive B reconstruction 
method. This method was originally developed for the SLD Rb measurement and is described 
in detail in 0. Briefly, the procedure is as follows. Well measured tracks with vertex detector 
hits are selected. In each hemisphere, secondary vertices are formed from these "quality" 
tracks using a topological vertexing technique ||. The most significant of these vertices that 
is significantly displaced from the IP is chosen as the "seed" . The B flight direction is then 
defined by the line joining the primary vertex and this secondary vertex. Additional tracks 
are attached to the vertex if they satisfy the following criteria. 

• the 3D closest approach to the flight direction is < 1mm 

• the distance along the flight direction to this point, L, is >0.5 mm 

• the L/D > 0.25, where D is the secondary vertex decay distance. 

These tracks, seed plus attached, are called "5-tracks" and a 6-tag is formed by calculat- 
ing their invariant mass assuming that each has ± ■ A correction is applied to account for 
neutrals and missing tracks which is based on the total transverse momentum of the vertex 
tracks to the B flight direction. Figure [l] shows a histogram of this "Pt corrected mass" , 
Mp T . Finally a cut on M Pt is applied at 2 GeV/c 2 . In Monte Carlo studies, this b selection 
method had an efficiency of 35% and a purity of 98% in '93-'95 and an efficiency of 50% and 
a purity of 98% in '96-'98. 

In addition to the Mp T cut, the set of "5-tracks" can also be used to study the structure 
of the B decay. This is done by fitting all 5-tracks to a single vertex and calculating the fit 
probability. Due to the finite charm lifetime, B decays with open charm will tend to have 
lower fit probability than those that are "charmless". Decays with fit probability greater 
than 0.05 are called "1-Vertex" and the remaining decays are called "2- Vertex". Table ^| 
shows that, indeed, decays without open charm are more likely to be in the 1- Vertex sample. 

A check is also performed on data in B — > J/^/X data where J/\I/ — > In these events, 

the 1-Vertex fraction is found to be 0.733 ± 0.094, confirming the efficiency estimated from 
the Monte Carlo. 
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P t Corrected Mass (VXD3 96) 

1000 r 




Mass (GeV/c 2 ) 



Figure 1: Pt corrected mass for '96 data and Monte Carlo. Cutting at Mp T > 2GeV/c 2 gives a 
very pure and efficient b tag. 

B Decay Mode 1-Vertex Fraction 

Single Charm 0.33 ± 0.01 

Double Charm 0.16 ± 0.02 

Charmonium + X 0.76 ± 0.03 

sg Model 0.73 ± 0.01 



Table 2: Monte Carlo estimates of the fraction of B decays of different types that satisfy 
the 1- Vertex cut. 

2 Inclusive Search for b^sg 

In the Standard Model, b^sg occurs through gluonic penguins and is expected to have a 
total branching ratio of approximately 1% ||. However, if the branching ratio were enhanced 
up to ~ 10 % by some non-Standard Model mechanism, it would nicely explain the "puzzles" 
described in section |l] Experimentally, it has been difficult to exclude even such a large 
B(b—>sg) due to the lack of a clean signature for these decays. 

The SLD analysis |10| searches for such a large enhancement by examining the high p t 
part of the spectrum, where pt is the momentum transverse to the B flight direction. 
Naively, one would expect K ±] s produced by b^sg to have a stiffer spectrum than those 
produced from standard B decays since they come from a direct b— »s transition rather than 



cascade b^c^s transitions. In a simple, jetset |TT] inspired model |L2[ it has been verified 
that this is indeed the case. The model, however, is quite sensitive to the choice of tuning of 
jetset. Table |^ shows the number of high-p t K s expected for different choices of tuning. 
"DELPHI Tuning" will be used for the rest of the analysis, but the tuning sensitivity means 
that limits on b—>sg can be set only within the context of a particular tuning choice. 

In addition to the M Pt cut described above, the 1- vertex cut is also used. Since b-^sg 
events are charmless, they are expected to be mostly one 1-vertex. Using the 1-vertex cut 
thus provides background rejection of standard b— >c decays as well as providing a 2- vertex 
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Tuning Choice 


K ± ,p t > 1.8 GeV/c per B 


jetset default 
DELPHI tuning 
No Parton showers 


6.9 x 10" 3 
10.6 x lO- 3 
13.8 x 10" 3 


b^c 


2.4 x 1(T 3 



Table 3: The number of high p t K s expected from b—^sg for a number of different choices 
of JETSET tunings. 



sample, which should not have much b—>sg in it and can be used to check the background 
calculation. The b^-sg signal should therefore show up as an enhancement in high p t K^s 
in the 1-Vertex sample. Only the high p t part of the spectrum is used both because it has 
high signal to background and also because the background in that region is well understood, 
leading to a small systematic error. Figure § shows the spectrum of K s observed in both 
the 1- and 2- Vertex samples for a sample of 50623 inclusively reconstructed 5's. Table f| 
shows the number of events expected from b— >c background and the number observed in 
the 1- Vertex sample and for all events. The data is well described by the & — ^ c Monte Carlo 
with only a small excess of high p t K ±, s observed. 



93-98 K 1 ,! Vertex 




0.5 1 1.5 2 2.5 1.5 2 2.5 

p T , measured Pj, measured 

93-98 ¥r,2 Vertex 




0.5 1 1.5 2 2.5 1.5 2 2.5 

p T , measured p T , measured 



Figure 2: The pt spectrum of identified K^s in the 1- and 2- Vertex samples. The data is 
compared to b^c MC background and to b^c background plus a 10% b^sg contribution. 
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K*^ > 1.8 GeV/c 


1 vertex 


All 




Raw 


Normalized 


Raw 


Normalized 


Data 


53.0 


1.46 x 10~ 3 


150.0 


4.18 x 10" 3 


b^c MC 


48.2 


1.29 x 10" 3 


132.3 


3.54 x 10~ 3 


Difference 


4.8 ± 7.3 


(1.7 ±2.6) x lCT 4 


17.7 ± 12.2 


(6.4 ±4.4) x 10" 4 


BR(&^s#) = 10% 


31.5 


1.13 x icr 3 


42.2 


1.52 x 10~ 3 



Table 4: The number K s with pt > 1.8 GeV/c observed in data and expected from the 
b^c Monte Carlo with and without the 1- Vertex cut. Normalized numbers are per B decay 
and are corrected for K identification efficiency 



The systematic error on the background rate is calculated by splitting the background 
high p t K ±J s into their component sources and assigning an error to each one of these sources. 
The largest source of background is from K ±, s from D 0, s that were produced in 5-decays. 
This source accounts for 23.5 of the expected 48.2 background events. Fortunately, the 
spectrum of such D 0, s has been rather well measured [13 and the SLD 5-decay model has 
been tuned to match it. This leads to a rather small systematic error of 4.3 events from this 
source. Other large sources of systematic error include uncertainties in the branching ratios 
of Cabibbo suppressed B decays of the type W^vs (3.6 events), p t smearing (4.5 events) 
and identification efficiency (3.3 events). 

The preliminary result is that the excess of 1-vertex K ±, s with p t > 1.8 GeV/c is 4.8 ± 
7.3(stat.) ± 8.9(syst.) events. Within the context of the b^-sg model with DELPHI tuning, 
this corresponds to B(b-^sg) = 0.015 ±0.023 ±0.028. Adding the statistical and systematic 
errors in quadrature, this gives B(b^sg) < 0.076 at 90% confidence. 



3 Measurement of B(B DDX) 

Measurement of the B "double-charm" branching ratio (B(B — > DDX)) can help to resolve 
the puzzles described section |l[ The SLD measurement of B(B — > DDX) uses the same 
inclusive reconstruction method described above to select a set of 5-tracks. The 2- Vertex 
cut is then applied - single vertex fit probability less than 0.05. The 5-tracks in these decays 
are then split into two vertices using a x 2 minimization procedure. The upstream vertex is 
called the "5" -vertex and the downstream vertex is called the ".D" -vertex. For single charm 
B decays, the assignment of tracks to these vertices will tend to be correct - the 5-vertex 
will contain mostly tracks coming directly from the B decay and the .D-vertex will contain 
mostly tracks coming from a subsequent D decay. However, for double charm B decays, the 
B-vertex will contain a mix of B and D tracks. Since come primarily from D-decays, 
if a K ± is found in the B- vertex, it is a good indication that the decay was in fact a double 
charm decay. 

To extract B(B — > DDX) from the data, the ratio Rk is formed: 

Rk = N K>B /N KjD (3) 
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N K - >B N K+tB N K - >D N K+tD N K . B /N KiD 

Data 1193 871 1830 1171 0.688± 0.020 

M.C. DD 3068 2303 2363 1865 1.270± 0.026 

M.C. tl Not-DD" 5301 3903 11387 6432 0.517± 0.007 

Table 5: Number of K ±, s in the B- and D-vertices for data and for Monte Carlo. 



where Nk,b is the number of K b observed in B-vertices and Nk,d is number observed in 
D-vertices. B(B — > DDX)may then be solved for by comparing R K observed in data to Rk 
found in Monte Carlo samples of pure double-charm and non-double-charm B decays. Table 
[5] shows the numbers that go into this calculation. 

Two corrections are applied to the data: background from light quark events is subtracted 
(AB(B — > DDX) = -0.9%) and the difference in efficiency of the 2- Vertex cut between one- 
charm and two-charm decays is corrected for (AB(B — > DDX)=-5.8%). 

The largest detector related systematic error comes from misassignment of tracks between 
B- and D-vertices. This effect is calibrated in the data using an initial state tag to identify 
the expected charge of leptons coming directly from the B decay. "Right" sign leptons 
should then be found in the 5-vertex and "wrong" sign leptons should be found in the D- 
vertex. The efficiency for correct track assignment can then be extracted from the ratio of 
correct to incorrect assignment. Similarly, B — > D* + X decays, with exclusive D* + decays 
are used to identify a set of tracks of definite origin, which are also used to measure the track 
misassignment efficiency. The statistical error of the track assignment efficiency is then used 
to calculate the systematic error due to this source, which is found to be 2.2%. 

The largest physics systematics are related to 5-decay modelling and to B(D — > K ± X). 
These sources lead to uncertainties of 2.1% and 1.8% respectively. The preliminary result is 
then B(B -» DDX)= (16.2 ± 1.9(stat.) ± A.2(syst.))% 



4 Measurement of r 5 +, r B o and t b +/t b q 

Measurement of exclusive B lifetimes provides important information about I?-hadron decay 
dynamics. Deviations from the naive spectator model are expected to be small and the ratio 
t b +/t b o is expected to differ from unity by only about 10% jnj. Large deviations from unity 



could, for example, indicate a larger than expected value for T(b— >cud) [ |T5|j . 

The SLD analysis [16] uses the same inclusive reconstruction method as described in 



section |1.2| . In order to improve the charge purity, tracks which failed the initial quality 
cuts, but which still likely originated from the 5-decay are also included. The total charge 
of the 5-tracks, Q is then calculated and the reconstructed B's are split into a neutral sample 
(Q = 0) and a charged sample (Q=±l,2,3). Figure [3] shows the distribution of Q for data 
and Monte Carlo. Monte Carlo studies show that the ratio between B + and B^ decays in 
the charged sample is 1.55 (1.72) for VXD2 (VXD3). Similarly, the ratio between B^ and 
B+ in the neutral sample is 1.96 (2.24) for VXD 2 (VXD3). fl 

^Charge conjugation is implied throughout. 
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-4-3-2-10 1 2 3 4 

Reconstructed Charge (VXD3) 



Figure 3: Reconstructed charge, Q, for '97-'98 data and Monte Carlo. 

Since the precision of the measurement depends heavily on the charge reconstruction 
purity, several methods are used to enhance it. In Monte Carlo studies, it was found that 
the charge reconstruction purity depended on the reconstructed M Pt , since decays that are 
missing some tracks tend to have lower Mp T . Therefore, to enhance the charge purity, events 
are weighted as a function of Mp T . For charged decays, the polarized forward-backward 
asymmetry can be used to tag the b or b flavor of the hemisphere. The opposite hemisphere 
jet charge also provides similar information. If the charge of the decay agrees with that 
expected from these tags, the decay is weighted more heavily. Conversely, if the charge 
disagrees, the decay is de-weighted. 

The B + and lifetimes are then extracted with a simultaneous binned maximum like- 
lihood fit to the decay length distributions for charged and neutral samples. Figure |] shows 
these decay length distributions for data taken in '97-'98. 

For the measurements of t b + and t b o, the dominant systematic error is related to b 
fragmentation. This is studied by varying both the mean fragmentation energy < x B > 
and the shape of the xe distribution. The systematic uncertainty from this source is found 
to be 0.036 ps for both tb+ and r B o. For the t b +/t b o measurement, the uncertainty due 
to fragmentation largely cancels out, since the two hadrons are assumed to have the same 
fragmentation function. For this measurement, then, the largest systematics are related to 
the fraction of 6-baryons produced (A(tb+/tb°) — 0.013) and B(B — > DDX) (A(r B + / r B o) = 
0.011). 

The combined '93-'98 preliminary results are then 

r B+ = 1.686 ± 0.025 ± 0.042 ps, 
r B o = 1.589 ± 0.026 ± 0.055 ps, 
t b+ /t b o = 1.061 ±° ° H ±0.027 

These measurements are the most statistically precise to date and confirm the expectation 
that the B + and B® lifetimes are nearly equal. 
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Decay Length (cm) Decay Length (cm) 

Figure 4: Decay length distributions for charged and neutral samples in '97-'98 data. 

5 Conclusion 

We have presented preliminary results on three analyses of B hadron decays. Within the con- 
text of the b^sg model described in [[12| we find B{b^sg) < 0.076, at 90% confidence. We 
have measured B(B —> DDX)= (16.2 ± 1.9 ± 4.2) %. And, we have measured the lifetimes 
t b+i t b° with the best statistical precision currently available. 

Work is continuing on each of the three analyses. The lifetime analysis will benefit from 
the addition of 150,000 more Z° decays from the Spring '98 running. All three analyses 
will benefit from a re-reconstruction of the data that will incorporate significant tracking 
improvements. And finally, more SLD running would significantly improve the errors of 
each of the three analyses. 
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